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<i>  n in  iv  *• » and  |iri'ii»fi!i'  response  w-ri-  aria  i yzcd  as  out-  problem. 
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in  tt  tiristi.  2 he  -e  ' i\J  i>.n  can  hi  rn'r-.i®-,  hut  tins  t . jit* 

v .» J 'ii.it  ion  would  -till  h<-,  i !i*i»er  s.isis  and  time  consuming, 
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i-  t !.i  vcond  ju.i  1 ■ . i s,  thi  pro  • e c l t I . ■ was  treated  as  .» 
r g id  • J.  , mil  it  s ;t.  i,-t  i it  |i.n  in  ,i  depth  >f  22  1 1 w.»  . calcu- 
i 1 1 e»i  . ue  in-race  Jvi  eli  i at  I.  n «»a  • V*  • c ' - in  the  -.oil  and 
IJ<‘  in  the  shall.  Iho  iinil  velocity  was  1981  ip>.  This 

i-  a t > ! ! ic  tent  approach,  hut  it  will  till  be  verv  time 

, ui  u.  . to  calculate  an  entire  penetration  event  in  this 

.sann-r.  fortunately,  hi  aase  of  the  quasi  -steady  tate  nature 
o;  ;•  net  rat  i .>ri  proct-s  within  a nominal  1 »-  homogeneous  geolo- 
c ic  lioi  , it  is  not  neci  sary  to  caiculati  the  entire  penetra- 
tiifi,  p- inf  t at  ion  dy-j.imiv-.  and  critical  force  loading 

c-ui.i  itaon  - :i  the  proie-.  tile  i.m  be  determined  with  reasonable 
KiuMi  . he  i.ialvning  ,mi  t tin  initiil  imbedment  and  those 
p.-rio.l-  nh'-ii  th-  pro  r*v  1 1 1 1 is  enteimg  different  geologic 
la  - i r •. 
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ihiN  report  describes  a i e-carch  program  to  mu^t  p;.itc 
the  tva-  ; 1 • 1 1 1 1 v .>!  the  application  of  a tw  o-d  imeiis  n>n.i  1 
f i n 1 1 e-d  i f f erence  code  for  the  prediction  of  projectile  pea 
etration  into  < irtn  ncdia.  It  was  conducted  for  th*j  s. 

\ * r i y ! hi  nicfi  Aat  erwuv:  eper  iment  Station  (K'ES)  *<y  Caliit;* 

no  hoi  to  ani  i ichr  ole\i  y , Inc.  f CR’i ) , under  Con*  rai  t HAC  t 
ja  . “4  - C • O'o;7  . ibis  program  was  sponsored  by  the  Defense 
Nuclear  Agency  ( I»V\  1 undt  r Subtask  SB20.),  Kork  Unit  la, 
"Contract  Support  Studies  of  i. round  Shock  Propagation 
tnrrugh  Soils  anJ  Rock". 

, 1 he  principal  n:\e-.tigutor-at  CRT  was  M.  H.  hagner  . 

I e<.  hn  i s a 1 advice  and  mdirue  wen  contributed  h\  k.  Ni . 

U;'i  nh.igen.  h.  S.  < ei  kc  perform  d the  computer  program 
- *ng  in..;  code  devt-loj  -tent  work  . t. . e.  I u*  ton  atd  (1.  !. . I.ggun 
; t v ’dec  valuable  ass i 't.mce  in  formulation  of  the  material 
models  and  data  reduction,  and  in  the  computer  production 
runs,  respectively. 

Tin  penetrator  spec i f icat ions , . s i te  profile,  and  mate- 
rial properties  were  provid'd  by  WES. 

The  Contracting  Uffi  o r 1 Representative  was  p.  V. 

Butler  of  hi s.  The  contract  ca  monitored  and  coordinated 
with  the  overall  UNA  penetration  program  under  DMA  Subtask 
sHdil,  Aork  'Jnit  PA,  "Contract  Support  Studies  of  Projectile 
Penetration  and  Impact".  \.  J.  Breithaupt  was  the  Contract 
ing  Officer  and  Col.  c..  U.  Hilt  was  the  Director  of  hi  S 
during  this  contract  investigation. 
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binite  .!  i f crer.ci  oar-ciival  mot  hods  h.»ve  the  potential 
for  overeor  i ‘ii..  the  i . m , t at  n>n  of  purely  malvtical  or  semi  - 
empirical  appra.i-he-  t >*r  the  following  rcaons: 

a.  The  dvnami.  mechanisms  of  target  response  and  of 
projectile  target  interaction-  can  be  explicitly 
incorpor  itc.l'  into  the  analysis. 


b.  The  basic  ; : per?  it*  of  target  media  (at  actual 

or  representative  sites)  can  be  incorporated  into 
the  analysis. 

C . I he  analysis  can  yield  time  historic"'  of  (1)  the 
■orce  J i,st  i i but  ion  applied' to  the  projectile  sur- 
face and  the  .oti'-equ  nt  projectile  deceleration, 
velocity,  and  position,  and  ( J ) the  stress  distri- 
bution and  defornat i >ns  within  tin  projectile  and 
its  internal  components. 


In  the  related  area  oi  ar  or  penetration,  finite- dif- 
ference numerical  technique-  h » e been  used,  beginning  in 
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I ht  ••  t • • : v\  • < i ti  v , ir  r«  i t r«  sea  rc  I have  tepri  to  ex- 
po t and  •••';»<!<  o;  t ho  pra  t 1 l .» 1 rob  1 cm  ■ s a apply  ing 

1 i>  finite-  i 1 1 • i > itv  • • 1 1 1 i an  - t o .i  ^ r.-.i.i  • • « - of  earth' 

pciM-t  t 1 1 ; >n  1 1 1 . 1*  m - , si...  * i-MMi  • t in  u t i t 1 1 v ■ f this 

!£•••:•  r »l  appro,!.  '«  % > i.j  i*  >i  d<  mo..  .t  rut  i*  n p roMcns  . 

!hi-  v*rk , whi.!*  wa  - primarily  arried  out  doring  the 

» i r halt  t 1 *4,  • t a i i ! s :.c  i -«  ;:h  ot  t ‘ ■ lomi  i.it  - mi  for 

the  oj»»  < , i- nt  pro.  ' m'  i in  a.’i.l  ; :ojrti!f  u pome  .m.ilr  e. 
**  ! i . h ,ir<-  d<  ■- . r i •>-.  d it:  it'..  I and  : . i»  o.u  b -w  1 , dge , 
for  « .«:»!-{  h , the  afiprua.  f,  ot  treating  the  pro_iec.:l»  as  a 
ii,  :J  bod  . ii.  t la-  fi  *CT«r  t e a 1 rtr.ai.'i  • of  earth  pom  t rat  ion 
d . : i • • . was  til'  u f i 1 i *ed • dur 1 ng  tills  [ r i a;i . H t is 
ir:  itt  it  lent  ft  t i r»!«ju*-  rhi.l,  r.d-  it  poi.  i It  to  first  ana- 
i>-e  ti.i-  vi  cr-.  * up*  t ii-tortton  associatiJ  with  penetra- 
tion '■>  rio-  pro  ie  til«,  a;  t then,  if  dr  i red , to  separately 
ami  vie  th<  re^pon-t*  ot  the  proto,  tile  itself. 


2 


APSK<hU'H 


i.: 


\ represent 
Me-  -a  i-n 

in ci. Ion.  e i:aj»  »c? 
till*  at  }';> 

so  i I o» , r *>  h i !<•. 


t*  hut 

hypothetical  i-.uti 

penef 

rat 

ion  pro 

a :».i  1 

- 1 - . ft  COf! ' (.ft-l 

o f t h e 

in 

mi! 

an  H 

. *1.  dot.  of.:  in  to 

e i .»g  i 

va  1 

pro  i ec 

o ,i  ; 

>*t  1 . sett  ’argft  1 

on  'id 

mg 

of 

\ a ■ 

• 1 St  lip,.  * 

AO  d 

in  t • 

KH*  s 

1 1 f i r. 

1 t e J ! 

5 t'or.-n,  e L a g r a n - 

g i an  . sit 

, V<  1 i i . , 

W .♦  ' 

1 

t '• 

the  h;i 

> 1 ■ to 

. i r m.i  t ec  hn  iqu<- . 

A.Ui-1.  IS 

hast’d 

the 

Hewp 

tr,<  t 

v.d.‘ 

i ViO 

types  of 

a tl.i  i 

V ' t ' 

of  t 

he  pr  ' 

! e c 1 1 ! 

■ pent  ? rat  ion  were 

pi  r f o rm--J 

! . In  th- 

t ! : 

>t  . t 

?!»?  p 

; ej ec  : 

ill  A , 

i . trtaic.:  as  an 

e l i st ; c ; 

■■'.i-f  , c 

I de 

i . > r-ri 

! c 

hud;.  . 

ihe 

tress  teponse  of 

tru  j • • • ec  tile  and  t : V J<  !■.  rr\i  * i on  o tin  r >.  rgc?  u«i  um  dur- 
ing j»eni.- 1 rat  ion  w».  re  contpwt  ej  a.  .«  ingle,  unit  icd  problem, 
fh  i s ~dl.it  ion  «.a  ,in  o .’  out  t ;i  the  miti.ii  ph  i i-  oh  pent** 
tr.it  u-n;  i.e.,  tint.!  t'e  pir-’t-c?  i It  no  e wa  s t ul  ly  embodied  . 
into  tin  earth. 

In  the  second  analyvi-,  the  so  i ! response  and  consequent 
force,  acting  on  tin  projectile  wore  computed  as  before  using 
the  regular  e i a t i c p.  I »>t  • c I Emulation  of  the  WAVE  -I.  code, 
hut  the  projectile  itseli  w.-m  assumed  ti'  hi  a rigid  body. 

In  typicil  earth  pem-T  ru’ ion  hruhleris  of  interest,  the  pro- 
jectile aiffi-r  oi  ly  eiu>-!i  trains,  and  these  would  not  he 
exreited  to  s i gn  i : ican’  ! \ inti  rac  t •»  i ? h the  ev«  . c 1 y - d i s tort  - 
mg  targe?  rt'd  i.«  ■ with  tin  ; onetrat  .on  ft  Cc  -os.  lor  pur- 
pM.es  (•;  ana!  't  mg  t he  j or.  tr.it. on  dynam  i r •»  and  t he  forces 
applied  to  the  pi  -sivti’e  surf;..  « , ft  i t h*-r  cfori  i -’asona  hie 
to  a .Mime  that  the  projt  tile  i a rigid  body.  This  approach 
h vp.a  >.sc-  the  detailed  c- 'input  at  i on  of  the  internal  pnxject  i lc 
ri'i-m  e,  there*  ;.  allowing  a rori  rap  ■ 1 and  efficient  calcu- 
i,i  t . <-n  of  the  ; - net  rut  i ->ri . Ti-:;  rigid  body  solution  covered 
penetration  through  tie-  i.  las.  r and  inti'  the  sahle,  to  a 
total  depth  of  12  If. 


f.f  '.-•mmakv  »>•  «t. -n;  is 

At;  t . . - . t .*»  pen.fr.it,  n of  the  e i a •«  f i . pi  a 'tit,  or 

de!  u~!  i‘  {■  : <ni . » reject  tie  w i > urricJ  out  for  »*.  s r»  >»*c  , or 
to  a Jer  t ’•  •■:  _ > 1 n,.  tie ■.  into  : In  >ml  layer.  The  n m i na  i 
celt  -i  in  tin  t.*rget  f>r  fhi  . saint  imi  was  I x I inch. 

Aff.r  p n.fr  »t  in,;  ib  at  (>  <n  f,e  , the  tieiage  dece lerat  ion 
r.  atii-U  aitoat  ss.,  g*-.,  ,bi  It  u l per- isted  tor  tin  balance 
c:  t nr  ,uf  :or,,  i n.I : . at  i ng  fnat  interactions  between  the 
target  oil  ana  the  tii-l  few  iiuhes  oi  the  nost  tip  essenti- 
ally determined  t fie  pi  net  rat  mu  dvriimi>_s.  The  maximum 
stresse-  in  tin.  sfei  j ;rcai  till  niurrid  near  the  nosetip, 
which  per  ten  ec,  iho.it  5.  a *.b,  .,rd  in  the  connecting  section 
'■v  *u*reti  *’■  n -•  and  the  *>odv  .•*  the  p rin  i c t i ! e , • uh  r c h exper  - 
i ••  "i  e.i  r -.til  c i K ‘ . 

bar  g * * cii!ati«ns  ucuin  i .1  in  the  force  history  applied 
to  t : >•  * «•  rr  a1- 1 • pr  ojei  ' i le  due  apparent  I y to  amplifying, 

cnchi  * mi'  in*’  i.ict.ori'  hetweeti  the  finite  nature  of  the 
co-pu'a* s rial  gild  and  t*>  -trmglv  liysf  eret  »c  load  unload 
b\  \ *•.  i r given  ny  f hi  soil  model.  ,\  method  for  alleviating 
thi-  problem  wa  found  during  the  rigid  body  solution.)  De- 
spite f!i  , oscillation  proi  Jem,  the  deformahle  body  penetra- 
tion .oiutieri  demon-.*  rate-  the  feast;  ility  ot  obtaining  a un- 
ited : of  ; s n a;  project  » It-  pern  t rat  : on  and  projectile  re* 
spon-e.  However,  th»-  i a vunheru'ic  and  time-consuming  ap- 
proach. A mort  practical  approach  i*  to  analyse  the  penetra- 
tion Jyr.ante'  and  projectile  respons  as  separi.  t«  problems. 
Iron  an  anaiy  i of  * he  penetration  of  a rigid  projectile, 
the  .surface  loaiing  on  the  projectil  - '•or  face  can  first  be 
determined.  let  loading  can  then  he  used  as  a boundary  con- 
dition for  analysis  of  the  response  of  the  elastic-plastic 
projectile  i t -»t  1 r . ihis  approach  ha  been  used  in  Ref.  5. 

The  rigid  body  penetration  solution  was  carried  out  for 
1 1 , ;>  ihmu  boring  this  time,  the  projectile  penetrated 
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rhrough  1 1>  -/'  :t  of  >oil  and  about  i ft  of  shale.  The  nom 
rr.il  cell  silt  in  t!u-  target  was  1x1  inihc-.'.,  The  average 
deceleration  level  wi  about  390  g’s  m the  soil  and  4 4 (»  g’s 
in  the  "iialc.  i he  discrepancy  between  deceleration'  m the 
soil  laver  ft  or:  t’he  .lefornahle  body  and  rigid-  bod  y penetra 
t : on  ■ - ' 1 at  : on  - probably  due  to  the  coarser  zoning  in  the 

rigid-  h •ii  v v a>c  . 

Targe  » * c.  r i l .i  t i ■ * n s i;,  the  force  history  applied  to  the 
pro!  c.  tile  a l so  occurred  *h  the  rigid-body  solution,  but  it 
wa  t:ut  t'hv-e  were  substantially  alleviated  by  using 

a :;oovot  i ; . a 1 i > "predef  ormed”  gr  id  in  the  oil  ahead  of  the 
project  lie.  this  technique  eliminated  ynchronour  interac- 
tion and  reduced  th>‘  «h  -tort  ion  of  cel’s  sliding  out  and 
ar  sind  * he  nose  tip.  In  trongly  hvsteri  t ic  target  media, 
•vohs  i 1c  rah ! c care  m.  • d - to  he  taken  to  realistically  model 
trie  uni  ruling  arid  reloading  behavior. 

f he  rig  id- body  penetration  solution  demonstrates  that 
f in l te - d i f f erenc e code  solutions  of  deep  penetrations  are 
possible,  although  the  cost  for  analyzing  the  complete  pene- 
trations,of  a high-velocity  projectile  into  a soft  earth  tar- 
get i>  I ikely  to  b...  high.  Fortunately,  it  t,s  not  necessary 
to  analyte  t he  complete  penetration  in  ordet  to  determine  ■ 
critical  design  condition  . In  deep  penetrations,  the  target 
resp-  mi  per.  '.ration  processes  generally  change  only 

gr  id  •,..»]  1/  with  irjcre.is  mg  depth  within  a homog  rreous  geologic 
layer.  :ipder  these  circumstances,  che  penetration  dynamics 
and  the  critical  venditions  of  force  loading  on  a projectile 
can  be  determined  with  reasonable  accuracy  by  analyzing  just 
the  initial  ri.--.~i  embedme r. . period  and  the  periods  where  the 
pt  oji-ctiie  no  e.>  is  entering  different  geologic  layers.  This 
appi  aac!.  ha-  been  used  in  Reference  A. 
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f;.  a!::,lit  ;»iis  tor  t h>-  pc  net  rat  ion- solut ion-  are  sum- 
c..;  r , c eu  tti  * ; v ; r . ! . 

2.  i 11  M ! k Vi  Oi 

f,t-  projectile  t ->  -bom  in  detail  in  figure  2.  It 

t >n  ; - -if  a :«»» I i < •».  c.  I : nd r i va 1 stxll  behind  an  ogival 

ri  Nu'tma.  I he  projc-tile  cavity  contains  a 2‘)5-lh 
-tool  I illa-t  < which  i-  connected  the  outer  shell 
i..  oiii..'ut  ric  i ni.,s.  ih-  fvtal  weight  i t>n2  lb,  and  the 
weight  lira  ; ft  \|  fitl  12  ps  i . 

i : t he  J i 1 1 ,i  b 1 1 ■ ’ o J v calculation,  the  steel  pro* 

net  i .v  f - treite.l'a-v  in  e i a - 1 i c perfectly  plastic  mater- 

ial hivin.;  tie.  propet  tic  - it. 'ted  in  Table  1 . The  equation 
of  state  and  other  detail-  of  the  model  for  steel  are  con- 
tained in  \ppe*u:  i . 

her  tiie  rigid  nod  v calculation,  the  only  information 
rtqiirc  l t ■ characterize  t he,  pro; ec t i 1 e was  the  specifica- 
tion of  fin  'liter  c»nt«<iir  and  the  total  weight. 

2,2  S I'll: 

The  hypothetical,  target  site  profile  consisted  of  a 
lt»  2/3.  *!  thief  liver  -i  -.oil  over  shale.  Set°  of  proper- 
ty curve  for  -e  media  were  provided  by  Waterways  Hxper 
irent  station  iW(S).’  '!  ia-se  curves  were  fit  by  the  hystcr 

e t i ; , i -t  it  - Ually  pi.  i-  tic  material  models  which  arc 
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described  in  Appendix  A.  Figures  3 and  4 show  show  the 
uniaxial  strain  load-unload  paths  provided  by  WFS,,  to- 
gether with  the  tits  given  by  the  models. 

The  models  are  seen  to  follow  the  WI:S  curves  quite 
well,  with  the  exception  of  the  f ina i • port  ion  of  unloading 
in  the  strongly  hysteretic  Layer  1.  Figure  5 shows  the 
shear  failure  surfaces.  (The  upper  plot  gives  the  low 
pressure  region  in  detail.) 

Table  1 summarizes  the  site  media  properties. 
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SECTION  III 


COMPUTATIONAL  METHOD 


3.1  HAVE  - L CODE  , 

The  WA\E-1.  code  was  employed  for  these  calculations, 
WAVL-L  is  a two  dimensional  code  which  solves  the  equations 
of  motion  tor  elast  ic -plastic  bodies  by  means  of  a finite* 
difference  Lagrang ian-ce 1 1 technique.  The  mathematical 
formulation  is  basically  the  same  as  that  described  by 
Wilkins."  A s, tiding  interface  was  used  to  allow  penetra; 
tion  of  the  projectile  through  the  computational  grid, 
for  these  calculations,  the  project i le/target  interface 
was  assumed  to  be  frictionless.  A lithostatic  field  in 
the  soil  was  generated  hv  prescribing  initially  compressed 
soil  states. 

5.’  GRID  DECOUPLING 

The  sliding  interface  formulation  in  KAVE-L  includes 
the  capability  for  dynamically  decoupling  or  unlocking  the 
grid  points  when  a prescribed  decoupling  or  failure  criter- 
ion is  met.  In  a penetration  problem,  decoupling  permits 
the  cells  initially  attached  to  the  axis  to  flow  around 
the  advancing  nose  shape,  sliding  along  the  interface.  In 
these  calculations,  decoupling  occurred  when  either  the 
generalised  plastic  strain,  , in  a cell  adjacent  to 

the  axis,  or  the  tensile  stress,  -o  , normal  to  the  axis 
(slide  lincj  exceeded  critical  values.  Generalized  plastic 
strain  is  defined*  as 


;i  ml  i ■»  .1  measure  uf  the  ^umiil.it  n't-  plastic  deformation  ;in 
element  undergo  o . d ■ ■ ' is  always  it  , and  thus  ^ 
mptiotonsoallv  i ut  reuses  as  any  type  ot'  plastic  strain  is 
experienced. 

■ f he  cittical  generalized  plast  u'  strain  and  tensile 
>tres  ■ values  u ed  tor  1.  «ver  1 and  Layer  1 were: 


Layer  I t -o  i 1 r—  .0?  — — 4 <)  psi 

Layer  2 (shale)  .013  12*'  psi 


Purine  1. 1*;;  - vei  oc  i ty  penetration,  the  soil  flowing  around 
the  nose  often  separate  iron  the  projectile  before  reaching 
the  tangency  point  ot  tne  nose  and  cylindrical  body.  i.. 
these  calculations,  separation  of  a point  could  occur  when 
the  stress  normal  tc  the  projectile  was  no  longer  compressive; 
i.o *'  i)  . When  thi-  criterion'  is  met,  the  point  is 
moved  as  a free  surface  point  and  may  subsequently  collide 
with  the  penetrator,  in  which  case  it  again  becomes  attached 
to  the  projectile.  (Closing  of  the  void  would  generally  be 
expected  during  the  final  stages  of  deep  penetration.) 

3.3  COMPUTATION  OF  STRI  - s(  s ACT  INC.  ON  PRO.IHCT  1 l.L 

The  projectile  contour  is  defined  by  a series  of 
boundary  points  spaced  along  its  periphery.  (In  an  elastic- 
plastic,  or  deformable  projectile,  these  boundary  points 
are  corners  of  the  cells  used  to  describe  the  projectile.) 

The  projectile  surface  lino  constitutes  a slide  line, 
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s t r e * e n art*  t ran-m  1 1 teJ  across  tie  slide  line  in  accor- 
dance i»ith  j friction  rule.  l-or  noth  of  the  present  cal 
cul  if  new,  the  surf ace  was  assumed  to  he  f r ict  imili-ss , so 
on  1 v ikimi.iI  stii'-es  wer*.  transmitted  across  the  slide 
line.  'he  normal  stress  applied  to  each  surface  segment 
of  the  projectile  is  determined  as  follows,*  Referring  to 
the  example  shown  in  the  --ketch  below,  the  stress  compo- 
nent-- in  i direction  norm'  to  t',.*  projectile  surface  ser 
•nent  - 5,  are  computed  at  the  centroids  a and  h of  the 
target  cells  whicn  are  in  contact  with  the  proiectile  sur- 
i -k  e segment.  , In  a i.agrangian  code st  res  ses  are  cell 
centered  quantities.)  By  linear  interpolation,  the  stress 


Pro  jt-c  t i lv 
surf  .ice 


Slide  line  between 
projectile  and  target 


Xbi.se  tor  of_^ 
surface  line 
segment  2-3  v 


■.Target!  cells  I and  II 


' Centroids  c">f  cells 


i 


f 

is  determined  at  point  c,  which  is  the  intersection  of  ah  S 

with  the  perpendicular  bisector  of  the  projectile  surface  ' 

segment  2-3.  I he  component  of  stress  at  c which  is  nor-  ’ 

mal  to  23  then  applied  all  along  2-3. 

An  alternative  method  would  he  to  apply  to  projectile 
surface  segment  J 3 the  ^t  res?  calculated  for  cell  I over 
segment  2 -it  and  the  stress  for  cell  11  over  segment  B-3. 

This  method  becomes  more  inaccurate,  however,  when  the  tar- 
get cells  become  sever  1 y d i storted . 
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AUCTION  IV 

Pl.Ni.T*;  \ i ION  Of-  U AS  I 1C  1’i.AS'I  1C  ( DCf-OKMABI.K)  PROJECTILE 

In  this  s lution,  tee  penetration  processes  in  the 
target  .-*>  * I and  frit-  response  of  the  projectile  were  simul- 
taneous I > anal'. . e.i  for  t:.c  interval  from  impact  until  the 
ne-c  was  full.  embedded. 

l.I  (lit  1 b‘  Pi  SION 

The  coriputat  irna  I grid  for  analyzing  the  elastic- 
piastic,  or  Uctomahle-bodv  penetration  is  shown  in  Figure 
o.  ( I h • • plot  - h 'w  • a crap  let  e c ross - sect  ion  of  the  a*  1 
svnnetru  problem.  The  actual  grid  contained  only  cells 
on  one  iJe  of  the  axis.)  The  nominal  lone  size  of  the 
projectile  was  .1  inch  in  the  radial  direction  (i.e.,  8 
cells  across  rh-  radius)  and  1.0  inch  in  the  vertical  di- 
rection. In  t(.e  body  of  the  projectile,  the  vertical  zone 
dimension  was  gradually  increased  to  a maximum  of  1.7  in. 
at  the  rear  ot  t lie  projectile.  There  were  about  300  cells 
in  t lie  pi  oj  ec t i 1 e . 

In  the  target  meditr»,  uniform  i x 1 in.  cells  were 
used  out  to  lb  in.  radio  , and  down  to  20- in.  depth.  Be- 
yond thi-.  region,  the  radial  and  vertical  cell  dimensions 
gradually  increased  in  a 21  geometric  series.  The  overall 
grid  extended  to  26- in.  radius  and  38- in,  depth,  using 
about  000  cells  for  the  target.  (Because  of  the  high  com- 
pressibility and  low  wave  velocity  in  soil  Layer  1,  the 
presence  of  the  boundaries  was  not  felt  by  the  projectile 
during  the  0. 8-msec  duration  of  this  calculation.) 
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The  Kline  runnim  along  the  pent*!  rat  or  contour  and 
along  the  a a i s of  the  target  war  epee  if  Led  as  a slide  line 
vs'  the  perutrator  impinges  into  the  target,  this  dual  line 
defines  the  perutrator  soil  interlace  and  allows  slippage 
and/or  separation  to  occur. 

The  lattice  point"  along  the  axis  were  initially 
coupled,  and  the  p.ojectilc-  tip ’and  the  corresponding  soil 
point  were  locked  together.  Ihus,  before  the  soil  point 
at  the  surface  could  nave,  nut  along  the  projectile  nose 
surface,  it  hid  to  decouple  by  having  the  assoc  iated  soil 
ceil  fail  according  To  the  decoupling  criterion  described 
in  section  Iii.  at  t iiernorc- , since  the  points  were  initi- 
ally locked  together,  the  initial  velocity  of  the  soil 
point  was  equal  to  the  projectile  velocitv.  This  produced 
an  unrea 1 i st ieal 1 e high'  initial  stress  near  the  impact 
point.  While  the  high  str'e.ss  decayed  rapidly  and  proha’  1 y 
did  not  significantly  at  feet  the  overall  solution,  it  u.t  - • 
subsequently  decided  tiiat  it  would  be  more  realistic  to 
have  the  n>  l and  soil  point  at  the  origin  initially  de- 
coupled and  to  assign. zero  initial  velocity  to  the  soil 
point.  • 1 h i s procedure  was  useJ  in  the  rigid-body  penetra- 
tion solution  described  in  Section  V. 

t.2  DEFORMABLE  PROTECT  1 1,1.  SOLUTION  RESULTS 

The  deformable - proj ect  i 1 e penetration  solution  was 
carried  out  until  the  nose  was  completely  buried,  with  the 
tip  at  about  20  inches  depth.  This  took  0.837  msec.  The 
computational  grid,  particle  velocity  field,  and  principal 
stress  field  at  this  time  are  shown  in  Figures  7 to  9. 
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In  the  I,. i g i a »g  i an  grid  1 1' igure  T)  t the  s > mb«*  1 "x"  in 
a ci*  I f indicate-  that  material  jn  that  tell  ha  s experienced 
plastic  flow  or  failure  at  sn^c  t :m*  during  the  solution. 

A volume  of  f !.*  •«  il  extending  outward  to  about  3 projec- 

tile radii  and  ih  ad  of  the  nose  tip  about  1 radius  has 
failed.  She  pr<*  i ,-c  t 1 1 e response  wi  - entirely  elastic,  ex- 
cept for  light  plastic  del  ut  nation  (T^  * .35^,)  in  the 

thin  annular  sict-i.  n connecting:  the  nose  with  the  body. 

A seri*-  of  t hi  fu  Id  plots  obtained  from  tills  solu* 
tun  a - the  pr>*ie-tile  entered  the  soil  is  given  in  Appen- 
dix u , 

1 he  dcnimi.'-  f thi  overall  projectile  are  summarized 
in  figure-  1«»  t . 1..’.  whu.ii  •.  >w  tint  histories  of  t he 

total  axial  foi  a.  plu.i  to  the  pr<  jecttle,  and  resulting 
average  deceleration  and  vel'vit-,  and  depth  of  penetra- 
tion. After  about  .3  - c,  cor  re -pond  inn  to  <•  inches  of 

penetration,  t!.<  in  rag.  iccelcr  at  ion  level  had  reached 
about  5*3  g's.  liii-'  level  persisted  until  the  end  of  the 
solution,  de  pite  the  fact  that  separation  between  the 
-oil  and  pro  u till  surf.tci  lid  not  occur  until  about  ;5b 
msec  when  the  n i e tip  wa-  at  about  11  inches  depth.  This 
indicates  that  axial  fn lies  applied  by  the  soil  to  the 
first  b inchc*-  «»t  the  n » tip  wen*  of  picdnminate  impor- 
tance t ) the  | lit  dynamics  of  the  problem. 

i he  large  o-c illation  which  are  apparent  in  the 
force  and  tcce  iei  at  ton  historic-  in  1-igun  1 <*  are  due  to 
interaction-  between  tie-  sunt-  nature  of  the  computational 
coning  and  r h<  -tror..j  ly  hy-teritic  character  of  the  unload- 
ing model  u.-ed  f.  tii-cr:  « fin  i.a.er  1 -oil,  This  problem 
will  be  discuss. d forth. r in  conni.tiun  with  the  rigid- 
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h«*»lv  solution  in  so'tion  VI. 

Figure  13  Mimnar i :e-  the  peak  compressive  stresses 
experienced  in  the  projectile  near  its  surlace  (solid 
curv«i  ind  i long  it  . ax  i (dashed  curve!.  The  peak  values 
are  preKil'l.  overstated  due  to  the  numerical  oscillations 
m the  surlace  loading  on  the  projectile,  bu t they  show 
the  location  of  critical  stress  regions  in  the  projectile, 
(he  nighv-t  stresses  are  experienced  in  th**  first  4 inches 
of  the  nose  tip  and  in  the  thin  annular  region  connecting 
tie  no<e  wit!»  tin  bo.lv.  Relatively  modest  compressive 
tress  lev*  l arv  s,.,  n in  the  wall  of  the  projectile, 
some  1 1 n s i uii  i s txperieiued  a<ar  th.e  axis  of  the  projec- 
tile; a-  is  shown  hx  the  vectors  having  leftward  compon- 
• nt  s in  Figure  . 

Figvites  lt  ana  i .«  show  mean  normal  stress  vs  distance 
profiles  just  i n i Je  the  projectile  surface  at  a sequence 
ot  times  after  impact,  i.e.  / : j j . )/3  in  the  row  of  steel 
cell'  along  the  periphery  of  the  projectile.  The  progres 
sion  of  increasing  stress  along  the  projectile  is  seen.  All 
of  the  profiles  after  .1  msec  show  a sharp  drop  about  5 in. 
hack  f.'"m  the  nose  tip.  At  later  times,  the  mean  normal 
stress  peaks  ;i,,n:i  i •»  the  narrow  connection  between  the  nose 
and  body  (near  2n  truhes).  Figure  in  shows  mean  normal 
stre:--*  vs  distance  profiles  along  the  projectile  axis.  Ten- 
sion occur-  at  the  leading  surface  of  the  slug.  The  decel- 
erating force  on  tlu  cylindrical  slug  is  being  applied  only 
around  its  periphery.  This  produces  tension  in  the  unsup- 
ported fiont  surfuci. 

Stte-s,  velocity  and  displacement  histories  were  ob- 
tained at  several  stations  in  the  soil  and  within  the  pro- 
jectile. These  are  available  to  interested  readers  in  un- 
published form  at  MS,’ 
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Axial  Distance  from  Nose  Tip,  in. 

Figure  14.  Early-Time  Profiles  of  Mean  Normal  Stress  Just  Inside  Projectile  Wall 
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RICH*  BODY  SOLUTION 

In  a Jeep  penetration,  the  target  material  undergoes 
large  distortions,  while  the  projectile  does  not  generally 
undergo  anv  significant  plastic  deformation.  To  analyze 
the  long  time  processes  of  Jeep  projectile  penetration, 
including  the  projectile  response,  with  a cor  > 1 e t e elastic- 
pLastic  code  solution  vould  require  an  impractical  amount 
of  computer  time.  This  situation  can  he  alleviated  by 
treating  the  projectile  as  a rigid  body,  while  analyzing 
the  soil  response  anJ  consequent  forces  on  the  projectile 
surface  using  the  regular  elast  ic • p last ic  formulation  of 
the  code.  The  rigid  body  assumption  is  justified  so  long 
as  the  projectile  deformat  ions  are  not  large  enough  to 
significantly  interact  with  the  target  response  or  pene- 
tration dynamics. 

The  advantage  of  this  approach  is  seen  tn  the  following 
considerations.  Computational  time  is  determined  by  the 
total  number  of  cells  and  by  the  numer  of  integration  cyc.es 
in  a solution  fi.e.,  the  total  number  of  cel  1 -cycles) . The 
number  of  cycles  is  proportional  to  the  reciprocal  of  the  time 
-.top,  whuh  in  turn  is  set  by  the  time  it  takes  the  fastest 
wave  (generally  a dilator  o>n  il  or  sound  wave)  to  cross  the 

* v 

smallest  cell  dimension;  ,t  1 (—  ) . Bv  describing 

max  c m l n 

the  projectile  as  a rigid  body,  the  computational  time  can 
he  reduced,  since: 

(i)  Larger  cells  can  be  used  (the  small  cells 

required  to  resolve  details  of  the  projectile  re- 
sponse are  no  longer  needed) 
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(7i  Longer  tine  s t ep  • ran  be  used  (both  because  nf 
• tin-  laigcr  cell?-,  and  because  of  the  lower  sound 
speed'  ummIIv  a lated  with  soft  earth  mater- 
ial- a ' contrasted  t.o  projectile  materials). 

I--  o:  tin'  approach  doe*-  not  preclude  analysis  of 
t : i interior  sv.-sp.in-e  ej  the  projectile.  from  t he  rigid- 
hod.  olution,  stress  loading  histories  imposed  on  the  pro 
jectiie  siri  ict  are  determined.  I'his  information  can  then 
?u  d to  "drive"  i detormable  body  analysis  of  the  stres 

*.»\e  re  ; -n  e m tie-  penetrator  itself.  Tins  coupled  tech 

r.  : le  ha  been  si-  -stull.-  enp  I!,yeJ  on  another  program.5 

; <:r:  :»  id  sif\ 

.he  con put  it  : >mi.i 1 gri  I for  the  rigid  body  solution  is 
.'ho.-i  ; i in. or.  i ",  I !u  >.  ntour  of  the  penetrator  per  i • 
p h v • - •»  i r e - o 1 v ed  with  a total  of  35  1 a 1 1 i c e points, 

spaced  at  m.  intervals.  the  computational  grid  in  the 

target  .-.i-  divided  into  two  principal  regions,  as  shown  in 
! igure  1“.  ! r> n the  problem  axis  ef  symmetry  out  to  a 

rad  ill-  of  jo  in.,  uniform  toning  of  2 x 2- in.  cells  was 
used.  be  vend  2>  in.  radio.',  the  radial  cone  width  was  in- 
cr-  ised  in  lu..  -.tip'.  Hie  gr  id  extended  to  a radius  of  7 
...  In  this  outer  region,  the  vertical  cone  width  Ac 
a 1 - • * in.cri.j-ed  (with-  increasing  depth)  in  1U»  steps  until 
it  reacht  i a maximum  value  f 10  *c0.  This  vlesign  for  the 
'out  -r  region  con -erves  computational  cells,  without  signi- 
ficantly it  let  ting  the  solution. 

A.  in  the  deformable  body  case,  the  K-line  running 
along  tii-.  penetrator  c»nt>  ur  and  along  the  axis  of  the  tar 
gt  r „a  -lult  line.  An- <fher  sliJe  line  was  used  along 
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DEPTH  FT 


•II*  m.  S5t"-4,  RICIO  MOt  RCNCTRRTtCft  (NT*  MU/OMU  RITC 


RROIUS  FT 


Fig  are  17.  Initial  ' onf  1 g«.rnt  ? on  frr  PI  pin-T'oHv  Solution 


■ thi  boundary  befwt ■••n  the  inner  and  outer  gi  id  rt'^iuii'i  in 
t«>  accommodate  the  ! i : *»er  ini;  a i i gnment  of  the  j*  i ids 
.*n  both  ' i d t* ' . ’«i  •»  « gn  i ! u ant  -lippage  occurs  at  This 

range.  This  type  -on  i iji»  has  !»eeu  a td  in  kefs.  5 and  4. 

A*-  di-'C.U'-i.d  in  s, t .<•!  ! . I , ? h<-  lattice  point'  at 
?hi  or. gin,  correspond  in.-,  t o t Ik  project  ill  tip  and  the 
sell  point  at  g round  tei-a,  were  nut  iallc  decoupled  for 
t f » • ricid  led.  '<*  lat  Hiit.  \ I t he  oil  points  are  thus 
ini r i . 1 1 v it  re s t . 

Kir.ia  i‘i<' 5 m.i' i m. i s •:  ;i  n»\  ri.s:?i  i\ 

i’ent  1 1 a * ion  i*  the  rigid  projectile  wa-  n eaer  u a I 1 \ 
ana  1 y *i-d  t r u..h  1 .« . i i i ' ll  and  for  ah>  ut  one  proiet  - 
tiie  length  into  i..r>*r  -ii.ilei,  to  a total  depth  ot  JJ..n 
ft.  i:u  gi  id  t .>u  i i y ,,  t . t i a:  it  an  i nt  c i *ued  i a t e penetration 
depth  d:  . 11.4  ft  a w;  it  : iguiu  1 h . I lie  "notches'*  in 

the  open  holt.  ht  h i.i.l  t t.t  projectile  are  due  t<>  pen.  die 
.deletions  of  d i - 1 • ■ » t » - 1 1 o.  11'  ifter  tin.  pi>>jectile  passes. 
Cells  in  this  i eg  i<*ii  no  longer  have  any  significant  effect 
on  the  penetr.it  ion- proce-  , and  their  rcreva  1 inprm  e* 
com.put  it  li'iial  eft  ic  h-ih  v. 

At  the  point  u.  the  .elution  shown  tn  figure  IS,  the 
wave  front  in  the  giound  had  nearly  reached  the  bottom  of 
thi  init,.«i  conput  at  i ori.s  1 grid  at  1 J . 2 St.  (Remember  that 
the  projectile  velocity  wis  supt  rsor.  it  relative  to  tie.  di* 
latational  wavi  vei»citv  in  tin  und  i - *uii  bed  soil.)  lo  ex- 
tend the  gild  without  inert  a- mg  the  nnnhi  r of  cells,  a 
"roll  ii'.!  grid"  piece  s wis  nH'd,  wherein  t tie  material  be- 
hind t'*-,  project  lie,  ii<'”>  the  ■air'iicf  down  to  5,7  ft  wa •. 
deleted,  and  matt  rial  ,»a  addt  d below  the  initial  grid. 
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Figure  2 i . Axial  Force  on  Penetrator  and  Deceleration  of  Projectile 
vs  Time,  Rigid  Body  Solution 
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23  and  24  at  early  time.'  and  at  times  after  about  ? msec 
(depths  below  about  13  ft). 

In  the  region  between  about  1 and  b.5  msec  (about 
2-12  ft  depth),  regular,  much  more  severe  oscillations 
at  about  b kHz  are  seen.  lhe.se  were  apparently  Jut  to 
synchronous,  amplilying  interactions  between  the  discrete 
nature  of  the  grid,  t tie  strongly  hysterctic  nature  of  the 
loading  tinloading- reload  tng  cycle  in  the  material  model, 
and  the  gross  distortions  of  the  cells  moving  up  the 
sides  of  the  projectile.  The  initial  loading  wave  velo- 
city in  the  Layer  1 model  about  540  fps.  Unloading 

(and  reloading)  always  occurred  at  a much  higher  vt  iocitv; 
i.e.,  of  the  order  of  bOOO- ?()()('  fps.  The  gross  distortion 
of  cells  moving  out  and  around  the  nose  caused  large  load 
ing  unload  i :ig  use  i 1 lot  ions  in  the  few  cells  which  wire 
applying  significant  axial  forces  to  the  nose.  The  pres- 
sure profiles  in  figure  25  illustrate  such  oscillations. 

lor  the  reasons  cited  above,  the  oscillations  in  the 
solution  between  about  1 msec  and  b . 5 msec  (2  ft  and  12  ft 
in  depth)  .ire  greatly  exaggerated  and  should  he  disregarded 
However,,  the  average  deceleration  in  this  region  ( 425  g’s) 
should  -till  he  fairly  realistic.  And  indeed  it  is  seen  in 
Figures  23  and  24  that  the  average  deceleration  level 'drops 
only  about  104  (to  about  3‘>0  g's)  when  the  severe  oscilla- 
tions go  away  with  the  introduction  of  the  prede formed,  or 
canted  grid.  Is  this  drop  due  to  the  grid  geometry  change 
as  ;uch,  or  to  the  elimination  of  exaggerated  forces  in 
the  synchronous  oscillations?  Apparently  the  latter:  the 
grid  canting  was  introduced  in  a very  gradual  fashion  ahead 
of  the  projectile,  and  it  required  only  the  very  slight 


change  in  grid  geometry  which  occurred  m the  first  few 
inches  *>e!ow  Id. 2 ft  (see  figure  1 ) to  eliminate  the 

synchronous  interactions.  As  soon  as  the  severe  oscil- 
lations were  reduced,  the  average  deceleration  level 
Iropped  to  5l>0  1 s . We  regard  this  as  the  more  accurate 

indication  ol  the  deceleration  level  in  l.aver  1. 


figure  dt>  give*  pressure  profiles  in  haver  1 (after 
introduction  of  the  predeformed  grid)  and  in  Layer  d,  and 
illustrates  the  smoother  prof i ! es  which  were  obtained. 

1 he  average  deceleration  increased  to  about  441)  g's 


in  the  haver  d shale.  I igures  17  and  d8  show  the  velocity 
and  depth  histories.  (A  constant  deceleration  ot  3‘.>U  g's 
was  assumed  between  d ft  irul  12  ft  depth.)  1 iguie  di1 
shows  the  velocity  vs  depth.  Inergy  partitioning  between 
the  projectile  and  target  dur  ing  the  course  of  this  solu- 
tion is  shown  mi  figure  3h. 


Stress, 
obtained  at 
These  are  a' 
form  at  WhS. 


velocity,  and  displacement  histories  were 
several  stations  within  the  target  media, 
a liable  to  interested  readers  in  unpublished 
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sin  ion  vii 

i'OSCl.US ION'S  ,\M>  kl-COMMliNDAT  IONS 

I he  results  nf  this  " t u J v provide  additional  informa 
t i*>n  with  which  • f **  evaluate  t lie  lole  of  finite-difference 
Ci>vie  ci  1 » u 1 a t 1 1 ■ u ' .H  earth  penetration  problems.  The 
f m i te- .1 1 f fereuce  approach  appears  to  he  very  useful  for 
anal.si-.  of  certain  ispi  cts  of  these  problems.  It  is  un- 
suitable,  or  at  lea  t ineflioient,  for  others. 


1 kl  .l!> 

ROih 

im;m 

1 KA  1 

IANS 

Sue ces 

-fui 

pern  t r a * < 

<r  deigns  will'  generally  exped- 

i  n c e 0 n . v e 

last 

i c d v 

: 0 im, 

itions.  (An  exception  might  be  a 

de  : gn  w i t h 

a -> 

aci  ! f 

1 c t a 

I ru-.o  section.)  It  is  therefore 

re  1 sonab  1 ■-  , 

and 

1 : f : 

cient,  to  anilyze  penetration  by  treat 

1 n t hi  , p r>* 

ject 

lie  a 

• J 

rigid  body.  Such  analyses  provide 

• nr ac  t i . a 1 

• . * 

h..J  : 

r : 

o ! camming  the  interact  ions  between  projectile 
shap-  , basic  tirgit  media  properties,  and  im- 
pact veli'cit;.  , (thereby,  providing  analytical 
guidance  fir  designing  and  evaluating  improved 
pro  ieet i 1 e des  igns I , and 

o |fet  e rm  i n i 1:4  the  external  loading  distributions 
on  the  projectile  (thereby  providing  bases 
lor  ana  1 vr mg  the  projectile  response  and  in- 
tern il  design  environment). 

Ihese  types  of  analyses  do  not  require  computer  solu- 
tions ct  complete  deep  penetration  problems.  They  require 
only  tbit  the  projectile  travel  far  enough  in  a given 

hi 


target  medium  to  reach  quas  i - steady- state  conditions. 

Such  analyses  should  therefore  be  yon  fined  to  those  por- 
tions o:  a penetration  event  where  the  loading  on  the  pro- 
jectile is  most  severe,  or  where  the  loading  changes  sig- 
nificant 1 >•  (e.g.,  during  initial  impact  and  embedment  of 
the  nose  into  the  surface  layer,  and  during  passage  be- 
tween different  geologic  layers). 

rhe  re  are  a nur.der  of  improvement  s which  should  be 
made  in  the  rigid  body  analysis  technique.  These  include 
efficient  tenoning  procedures  to  correct  grid  distortion 
in. the  target  medium  near  the  penetrator  surface,  and  the 

establi  -hment  of  minimum  coning  criteria  to  ensure  reason- 

ab!v  accurate  solutions.  Both  of  these  steps  are  being 
undertaken  in  a current  program. 

Because  of  the  iong  computing  times  required,,  it  is 
probable  inefficient  to  use  a f i ni ? e - d i f fe fence  code  to 
analyr*  complete  Jeep  penetration  problems  (i.e.,  from 
impact  until  ttu-  ,,ro;ectilc  comes  to  rest). 

7.  - , PR0JTC7  I Lh  KliSl’ONSi: 

fo  analyze  the  response  of  the  projectile,  a solution 
of  the  penetration  may  be  performed  in  which  the  projectile 
is  treated  as  a deformable- body  (such  as  was  done  in  this 
study).  However,  f i n i te  u i t I \rence  solutions  of  this  type, 
in  which  the  projectile  deformation  and  the  target  deform- 
ation are  treated  are  time  consuming  due  to  the 

markedly  different  character  of  the  projectile  and  target 
responses.  (iood  spatial  resolution  is  needed  in  the  pro- 
jectile in  order  to  resolve  structural  details,  and  the 
wave  speed s in  typical  projectile  mateiials  are  relatively 
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AIMM.M*  1 \ A 


MA  I I U I \1.  Mt'DI.l.S 

Material  rodeis  tor  the  two  site  layers  were  devel 
Oped,  based  o:  ti.e  set  >•:  iet.ro  <nt.it  tee  s,.(!  p rope  I t i c*  - 

tarnished  by  Wt  S.  Hvsteitic  elastic  ideally  plastic 
model-  u-itin  a 1‘randtl  Kiu'-s  i non  assoc  i at  < d I tlow  rule 
were  employed. 

\ non  hy.-teret  k e 1 a - 1 i c • pc r t ec t 1 v ■ pl.stic  model  was 
used  tor  the  steel  p: . i e<.  ' i 1 <•  . 


\.  ! 1 AVI.!'.  1 

A.l.l.  halation  d1^_  - tan 


1 he  pres  sure  i s g i ecu  by  a ‘lie  tlriinei  en  equation: 


where 


1*  - I a ..  j ♦ , e 

e i-  spec i tie  internal  energy 
h - pressure 
, density 


. - < I r Cine  i sen  constant 

lii  h yd  i'os  tat,  i . -pec  il  led  through  a bvsteret  Jc  model 


A. 1.2  Loading  Hydio  tat 

The  leading  h v <!  i t a t up  to  the  point  of  void  closure, 
u , was  tit  from  i’  ..  data  deduced  from  the  loading 
curve  of  stre-  vs  train  ami  stress  difference  vs 


/ 


pres. are.  Ik  v>rui  imJ  closure,  a we  1 1 - hcha  ved  asymptotic 


fe • I - enp  1 in  i d . 


* u--0 


\ . * h i'vjh  . ! 1 


. - 1 . * h 

V 

max 

< . 

) 

i K,„ 

■ * 

*'l  e\:> 

X 

■k  s 

V , i ,e  ■ f t 

Ml* 

^ • ‘I 

. " t .iu  t .i  ro  : 

\ - J .<_••!  1 

A 

U' 

1 Mli  « 

m a x 

. (nH'.Uo  Mb 

it  = 

A 

i * 1 

W-  « 

*ie;  ,z~ 

1 - J . II 1 

v 

1«>. 

1 •if  | r - 

l . i 

i ; 

. (;4i< ! A 1 M" 


. (>5<i4 


1 . '.I J s 2 gm/cm' 
{ » 2l>  Jh/ft3) 


shear  Moduiu-  _la>a_d_(j£ 

I lu-  s hea r ni  jul  i • function  fur  loading  was  determined 
i; , jpliical  a i sm  Tit  ot  t he  variation  of  Poisson's 
ratio  along  the  urio.iiny  -t  path  (stress  difference  vs 

prc-;,uri  I and  tin  'ulj.  n<kilu  (a>  derived  from  the  hydro 
,fa'  equation-].  \ : ,m  n:.  a I *n  equal  to  the  unloading 
spear  rndnla  was  rq  "u-.i  t<  pr-vitit  energy  generation. 


l’load 


38(1- 2 ■>) 

7 ( 1 * v)  ’ uunload 


B « fp*i)  d,Joad  i \ 

dli  i 1 ( 
> I 

v * ’5  •'  -2S  CXP 


Initial  values  of  the  bulk  modulus,  the  shear  modulus, 
the  constrained  modulus,  Poisson's  ratio,  and  the  dilata- 
tioral  sound  speed  from  the  model  are: 


Bo  3 4»“37  P*>  vQ  = .2 S 

Go  = 2-542  P*1  cQ  = 542.6  ft/sec 

Mq  s 7,t)2fc  psi 


A . 1 . 4 Unloading  Mod e 1 

In  the  hysteretic  region,  u< ac  , a model  giving  a 

fan  of  constant  slope  l’-u  unloading  paths  was  used 

(slopes  are  a function  of  u }• 

max  J 


F 


F 

max 


♦ K ( M - u ) 
v max7 


u <u 
H«ax  Mc 


= Fc  * *cU-uc) 


u >u 

max-  c 


The  unloading  hulk  modulus  varied  from  a minimum  of  100  ksi 
to  a maximum  of  702.8  ksi.  (The  maximum  corresponds  to  the 
bulk  mod u 1 u<  tor  loading  at  void  closure.)  The  unloading  bulk 
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modulus  «a>  made  a function  of  t he*  maximum  compression, 

“max  * ' constant  value  tor  the  unloading  shear  modulus 

wa  assumed.  li.e  initial  I’m  i sson  * s ratio  lor  unloading 

v‘;*'  thc  ",mc  us  tVr  I’^d.ng,  v - .25  , giving  0 

' SiSl'  I**  lir  t ( s the  above  equations,  \'  is 

the  hvJiostat  loading  p ft  oil re  corresponding  to  g and 

,•  . max 

fc  i s a linut  mil  given  nv : 


h,  * I-.  eip  (b„  ) 

. I ‘ “max' 


u <u 
max  Mc 


Ihe  constants,  as  determined  from  the  given  unload  paths, 
are: 

i'l  3 t» . H«>-I  r>  x H»  Mb 

b j - I'.'Jlii.  x 1 (l  Mb 

t = 2 1 8 . 2 


(or  ..  vc  . unloading  is  assumed  to  be  reversible 

i . e . , i * I o 2 ! i"a.  the  I o.id  path). 

os  tat  ii  tension  is  limited  by  imposing  a minimum 
value  ot  pressure:  I*  . ^ }(.  nsi 

min  • * • 


t . 1 . 5 l a 1 1 tire  S u rj  a e e 

A yield  function  V(l’)  ws  fit  to  the  failure  enve- 
lope supplied,  where  V \ at  failure.  The  functions 
derived  are: 


T *■  3,4473x10^  + . 4851.41'  - 2 . 3 5x  105P2 


I I , MIxl' 


l.""xl06(4.8droX10  f>  • Pr*1’15, 

0 - I'  • 4 . 8 2(>.>x Id  -'’ 

V = ’ 4 . 8 Jt>5x  1 »'  ■'  • I’  ■ 1 .3780*10’  3 

> - l.Jlixl'i  r'  ♦ .(»!*»’  3 (|>  1 . S'MDx  1 n ’ 'j 

i’  1 . 7 H f > 

Here  the  jrii,t  tor  N and  P are  Mh.  PI  >t>  of  : ht  failure 
surt.K'c  ate  i tie  » udetiv.  i n t tie  it- 1 low  in*;  set  ->f  »'  i gures . 

A.  1 . " Crrjp  i ra  t t v_e_  ke - id  t - 

Pint"  e t stress  ^ train.  s t re  - s • d i f to  react*  vs  pres- 
>ure,  ar<J  pre  -ate  t,  - c.> ! ame  t i i i strain  for  uniaxial  strain 
load-unload  p..t!.s  c * »»:> puted  with  the  modi),  along  with  the 
W!  S curve  , a.r  e hown  in  I-  igui  e 3 in  the-  text,  and  in  l:  i g - 
ares  A - J t \ ' ; rs  1 1. 1 ‘ a upend  i x . Plot.'  of  the  Hugoniot 
and  hydro  t.if  nj  ’ high  pr<  - -tires  ( SO  kb  t aie  shown  in 
I A- 4. 

A.:  LAY! H : 

Mu*  to  t lie  anticipated  lew  stress  env  i ronment  , . the 

pressure-  t >>i  1 «yer  d were  t a! en  direct  1 v from  the  hydro- 

staf,  w . t :i  no  energy  eon  t r i 'uif  ion  . 

A. 2.1  Loading  ‘-lodel 

The  following  equation-  were  deduced  for  the  loading 
h yd  r • ’ -tat t 


~ 1 


l>  * a; 


* * b 

w * U C 


= P-  * KB.4*‘  ' ‘lV..*U‘-*  I t* 


I*  , * ►* 

V, 


ITIJX  c 


•xp  — r 


U^O 

^<W2.MC 

W u „ 


lhe  values  of  The  constant.',  jre; 


a = 5.7348x10  *'  Mb 
b * 1 .3“i»»r,xl0  s Mb 


c * ■*  Mb 


.=  .OOolhS 
c 


K 

ma  \ 
»’ 


.048263  Mb 

.089465  Mb 

5 . 7 604x1 0 Mb 

2.0505  gm/cm3 
(128  lb/ ft3) 


The  shear  modulus  for  loading  was  specified  to  be  con- 
stant at  a value  of  20  ksi. 

In.tial  values  ef  the  hulk  modulus,  the  constrained 
modulus,  f’o  is  son’s  ratio,  and  the  dilatatiunal  sound  speed 
from  the  model  are: 

B = s 3.18  ksi  . v = .389 

o o 

M = lo;».8  ksi  . c = 1 ,994  ft/sec 

o o 

A . 2 . 2 IJn  leading  Model 

In  the  hysteretic  region,  a , a model  giving  a fan 

of  constant  slope  !’  . unloading  paths  was  used  (slopes 

are  a function  of  a >: 

uax 


where  ! is  obtained  Iron  the  loading  hvdrostat  ectua- 

;ua  v , 

t ion  ' evaluated  at  ..  ..  lac  constants,  as  determined 

III  d X 

from  the  niven  unload  path" , are: 

K = .01378'.'  Mb 

v 

Kj  --  . ■ 1 4 S J t> 3 Mb  • 

foi  a , unl  oidinc  is  assumed  to  be  reversible 
; 1 t : i ;.••.»  < r ho  ! oad  j.  i T i : . 

Hydrostatic  tension  is  limited  by  impo^-inn  a minimum 
value  of  pre:  sure : I'  = l'U  psi. 

A . 2 . 3 I a 1 1 tire  Sur  lhn  e 

A yield  tunctioi.  t(1’j  has  fit  to  the  failure  envelope 
suj; nested,  ..beta  V * .31 ' at  failure.  The  functions 
de  r ; v ed  a t e : 


r 


'i  s - 3. 527oxi n4 r 


-6.89  if»x  1 o'  *’•  !’*  0 


\ - . i “ ix  ii> 


't*.\  1 ir*  ( 4 , b Jb3x  1 0*  ^ - 1* ) *'  ’ l'’1’ 


0*  I*-- 4. 8:63x10 


> - . I 7 1 x 1 0 


4.8203xl05* p<j. 4131x10* J 


> »’ 3;  I?l.xl»"9  * . <*  1 53S  5(1’  2. 4131x10'  ')  I’-:.  4131x10 


-3 


lit  re  tht  unit-  tor  t and  I’  are  Mh. 


\ . 4 l »;:)]). i rat  it  e kt  -nit  s 

. l'l)t-  i : <r  s vs  -tram,  stres>  Jifterence  vs  pres- 
suu  , an  I pressure  volumetric  strain  for  uniaxial  strain 
load-unloaJ  path.-  computed  with  the  model,  along  with  the 
wl.S  curves,  are  shown,  in  I igure  4 in  the  text,  and  in  Figures 
A 3 arul  \-f»  m this  appendix. 

A . ~ STHd. 

\n  - ex:  d i ;;;4  equal  ior.  of  state  for  iron  was  employed.* 

' Appropriate  riexhas:  ncal  and  strength  properties,  for  the 
•-t'H'i  all  oe  ot  t h*  p.-net  rator  were  tlien  added  to  the  model. 


*M.  1?.  ■.vifjnt-r.  A-.  -• » ' ’-a  -U-e  i-r  • ' i l-  l •-j-'jk  a».d 

- . ..  • , she  k Hydrodynamics  Report  "384  <H> , August 

1 "1 
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The  following  properties  were  assumed: 


Normal  density,  , : 


Hulk  modulus,  K : 
o 


'.84  gm/cm3  (48-)  lb/ft3) 

! Mb  I J. 445x10  psi 


Shear  modulus,  .85  Mb  t 1 . dl>4.\  1 0 ps  i ) 

Yield  strength,  V:  .Ull544"5  Mb  ( SO , m»n  j.si) 


The  shear  modulus  was  assumed  to  he  constant.  The  bulk 

mod  n Si'  varies  through  tie  e.ju  »t  *«»n  of  state.  A conventional 
v ' n V t • ' . i e l d ’ ! w.:  - I n rj  ; j , e . , - t a t « s are  elastic 

: ' i A-*  *,  > c-'tis  t ant  . 1 1 J i o • t 1 1 : c tin.  ioj.  was  1 i m i ted  by 

: -![!'■'  i n<  ,i  v a 1 «e.  >>t  u t are: 


, -Is 

r:  i n 


hit!)  t fi : : condition,  the  tensile  .• : •*  cc  is  limited  to  a 
value  i . ’ 

I'ht  low -pressure  ( sub  t runs  i t ion)  portion  of  the  pre- 
v i v a - 1 . :■  t mu  I a ted  tu.ua  t ion  *f  state  of  iron  covers  the 
ra t\>'  o{  interest  and  is  given  !>v: 

I’  - J..  ( 1 e _<0 


i j “ ♦ a ..  I i 1 - . e 


0<ls- j 


1 he  value-,  of  the  constants  are: 


= I.nSn  Mb 


i = 1.81 


. J HSU  Mb 


= .IUl 


a.  - ?o.551  Mb 

a 


* 


inure  \ . . Uniaxial  Strain  Unloauinn  Paths  ami  I'ailuro  Surface  in  Hii;h -Pressure  hi  n ion 


H'kujv  Uniaxial  Str.nn  Strf<>  for  lave 
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IMSI'ORTIO.N,  lILiH-m,  AM)  Pk  1 NCI  PA  I STRESS 
I I PL!)  1*1  O 1 > D!  lOl.MABI.i  - BODY  SOLUTION 

A series  of  t ield  plots  of  the  comput.it  ional  grid, 
particle  velocities,  and  principal  stre-se-  which  illus- 
trate t fie  projectile  and  -.oil  response  during  the  de  lot  nah  1 e 
hedv  .-ol.it  u n welt  prep.it -.d.  Tlu  si  are  -nu,.n  fur  the  follow 
j Hi,  penetration  depth.-  ana  times. 


f-  Depth 

I me 

figure  No. 

of  Plot  for: 

(in.) 

! ,.  ec  ! 

Hi  stort  ion 
; ( ’>  r i d 1 

Ye  1 oc 1 1 v 

St  res s 

• r>(> 

d 5 . 3 
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The  velocity  vector  field  plots  shew  the  direction  and 
magnitude  of  the  particle  velocity  at  each  lattice  point 
in  the  computing  grid. 

I n tiie  st  res-  field  plots,  the  principal  components  of 

t t ! - * « ! ' » • s > 1 u • he  e -unveil!  {,.*i  -. ; 


TIu*  magnitude  of  the  two  principal  strokes  in  the  r-z 
plane  are  plotted  in  -flair  corresponding  principal  direc- 
tions. !he  thi'd  principal  stress  (in  the  azimuthal  di- 
rection) 1 >»  plotted  along  the  line  bisecting  the  other  two 
principal  directions.  Vectors'  pointing  to  the  right  are 
compressive,  to  the  let'*,  tensile.  An  example  of  how  a 
stress  tensor  is, plotted  is  sketched  below: 


The  first  three  figures  of  the  grid  (B-l,  B-4,  and 
3-5)  shew  the  initial  entry  of  the  penetrator  into  the  top 
surface  of  the  soil,  \ local  rczone  is  required,  as  shown 
in  figure  H (> , when  the  surface  grid  point  contacts  the 
peiiet  rat  or  surface  (see  figure  H-5). 

In  figure  B-7,  the  stress  loading  of  the  projectile 
nose  after  s"  penetration  is  seen.  Note  the  strong  verti- 
cal orientation  of  the  stress  vectors,  due  to  the  free  sur 
face  Jt  the  periphery. 

figure  H-1I  show  another  rezonc  which  was  performed 
to  rt.  orient  the  grid.  A small  portion  of  tite  distended 
soil  that  was  living  off  was  also  removed  at  this  time. 

figures  B ID  and  B-14  show  a buildup  of  radial  and 
•,oop  tension  along  the  iower  end  of  the  slug  and  a high 
vertical  compressive  stress  in  the  thin  section  connecting 
the  nose  and  the  slug. 

I he  detachment  of  the  soil  from  the  penetrator  and 
hole  format  ion  .may  he  seen  in  figures  B- 17  and  B-18. 
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APPENDIX  C 


DISTORTION,  VELOCITY,  AND  PRINCIPAL  STRESS 
MELD  PLOTS  - RIG  ID- BODY  SOLUTION 


A series  of  field  plots  of  the  computational  grid, 
particle  velocities,  and  principal  stresses  which  illus- 
trate the  soil  response  during  the  rigid-hodv  solution 
were  prepared.  These  are  shown  for  the  following  penetra- 
tion depth  and  times: 


Figure  No 

. ot  Plot 

for: 

Depth 
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(msec) 

Distortion 

(Grid) 
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Stress 

.34 

.17 
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G - 21 

22.31 

11.58 

C-  22 

G-  23 

G - 24 

The  plots  of  the  grid  also  show  which  material  is  cur- 
rently failing;  i.e.,  is  on  the  yield  surface.  This  is  de- 
noted by  cells  containing  an  x'  or  + , with  x indica- 

ting a compressive  pressure,  P >_  0 , and  + indicating 
hydrostatic  tension,  P < 0. 

The  velocity  vector  field  plots  show  the  direction  and 
magnitude  of  the  particle  velocity  at  each  lattice  point  in 


the  computing  grid.  For  clarity  in  viewing  the  soil .re- 
sponse, the  velocity  vectors  of  the  projectile  are  not 
shown. 

The  plotting  convention  for  the  principal  stress 
field  plots  is  the  same  as  that  used  in  the  deformabld-hody 
solution  (Appendix  B) . • • 
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Figure  < -14.  Particle  Velocity  Field  with  Projectile  .it  10  ft.  Depth 
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('-17.  Particle  Velocity  Field  with  Projectile  at  15  ft  Depth 
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Figure  (.-1°.  < rid  (.mf  lgur.it  ion  with  ('ro)c,  tile  at  16.9  ft  Depth 
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Figure  C-20.  Particle  Velocity  Field  with  Proloctite  at  16.9  ft  Depth 
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